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The reweighting scheme developed in Glasgow to circumvent the lattice action becoming complex at finite 
density suffers from a pathological onset transition thought to be due to the reweighting. We present a new 
reweighting scheme based on this approach in which we combine ensembles to alleviate the sampling bias we 
identify in the polynomial coefficients of the fugacity expansion. 



1. Introduction 

Lattice simulations of dense matter are impor- 
tant for the understanding of future heavy ion 
collider experiments as well as the equation of 
state of neutron stars [Q |2| . The main obstacle to 
the numerical evaluation of QCD in this regime is 
that the lattice action becomes complex with the 
inclusion of the chemical potential /i, prohibiting 
naive probabilistic Monte Carlo methods. Recent 
attempts to evaluate QCD at finite baryon den- 
sity have involved evaluating canonical ensembles 
with a given number of background quark sources 
(HQ, and the evaluation of QCD-like models in 
which the lattice action is real at finite /i 
An earlier reweighting approach was developed at 
Glasgow, in which the Grand Canonical Partition 
function is evaluated semi-analytically. 

2. The Glasgow Method 

The fi dependence of the lattice action is 
made analytic in the Glasgow reweighting method 
through the formulation of a characteristic poly- 
nomial in the fugacity variable z — exp(/i/T) 
. The fermion matrix M (defined with Kogut- 
Susskind fermions Q) is re-expressed in terms of 
the matrices which contain only links between lat- 
tice sites in the spatial directions G, and forward 
and backward in the time direction V and V'. 
which allows the definition of the propagator ma- 
trix P [§. 

2iM = 2im + G + Ve M + U t e" /i (1) 



The propagator matrix is then used to re- 
express detA/ as a characteristic polynomial in 
the variable e~ M , 

detM = det(G + 2im + + W) (3) 

= e n - n * ntf * det(P - e - ") (4) 

2n c n^nt 

= e n ° n « nt >* c n^ n11 (5) 

where n^nt is the lattice volume, n c the number of 
colours, and the expansion coefficients c„, so de- 
fined, are functionals of the lattice gauge fields. 
Since V is an overall factor of P the expansion 
can be further simplified through the symmetry 
Z, H associated with performing a unitary trans- 
formation on P by multiplying the timelinks by 
e 2mj/n t ^ wnere j j s an integer. This then allows 
relation of the expansion coefficients of the char- 
acteristic polynomial (where n t = 1/T) to the 
canonical partition functions Z n and the Grand 
Canonical Partition function Z(fi). As Lee and 
Yang showed with an fsing ferromagnetic system, 
in the thermodynamic limit a phase transition oc- 
curs wherever a zero of the fugacity polynomial 
approaches the real axis in the complex- z plane 
JTcj ] . The zeros a n are determined numerically by 
rootfinding the reweighted polynomial expansion 
coefficients. 

Z(fi) = J DU detM(fi) e~ s " (6) 

= z » e "" /T ( ? ) 
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\detAf( Mo )/ po (9) 

^ g-ncnSnt/i Jj ( e »*M _ q^) (10) 
n=l 

Reweighting introduces normalisation by the 
ensemble generated at into Eqn.(^), which 
leaves the analytic determination of the critical 
points unaffected. The advantage of the scheme 
is that, even though the fermionic action is in 
general complex in SU(3), an ensemble can be 
generated at \x Q = 0. However, it has been shown 
that the unphysical onset transition at /i = ^m^ 
of quenched measurements [|ll], persists despite 
the inclusion of dynamical quarks. From which 
it is concluded that the Monte Carlo sampling 
is ineffective, because the fi = ensemble has 
little overlap with the physically relevant region 
|l2j [^3) . Parallels can also be drawn between this 
pathology and the sign problem of the related 
reweighting of the Hubbard model |l4|] . 

3. SU(2) at Intermediate Coupling 

3.1. Composite Reweighting 

For SU(2) we can vary the value of /i we use to 
generate ensembles (as the group is pseudoreal for 
quarks in the fundamental representation p5[). 
This then allowed us to investigate the reliability 
of the ensemble-averaging of the polynomial ex- 
pansion coefficients of Eqn.(||), which is defined 
through the fugacity expansion. 
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These coefficients are evaluated recursively 
from the propagator matrix P, and the denomi- 
nator detM(/i G ) by evaluating Eqn.(|^) at /1 = /i Q . 
Where the ratio on the lefthand side of Eqn.(|l2|) 
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Table 1 

Dependence of r\\ for SU(2) at intermediate cou- 
pling on the fi (upper), and number of ensembles 
included in the composite reweighting (lower). 



differs from one, we found from our measurements 
that the reliability of the ensemble-averaging is 
strongly affected. In turn this affected our deter- 
mination of the Lee Yang zeros through rootfind- 
ing, and our measurement of related thermody- 
namic observables evaluated from the expansion 
coefficients. We monitored the effect by measur- 
ing the ratio W 21 (n), of the expansion coefficients 
from ensembles generated at two different values 

of [l Q . 
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Although this ratio should be independent of 
n, our measurements indicate that ensemble- 
averaging selectively affects the determination of 
the expansion. Only a small range of n can be 
reliably determined from an ensemble generated 
at a given value of \x . However, since this small 
range centers on different n as \x Q is varied, coeffi- 
cients from a covering of several ensembles (gen- 
erated at different values of /x ) can be combined 
to alleviate the bias inherent in the reweighting. 



The reliably ensemble-averaged range of coeffi- 
cients from the ensemble generated at \ii (multi- 
plied by W 21 ) is used to replace the corresponding 
range of the expansion of the ensemble generated 
at fi\. As the number of covering ensembles is 
increased these new composite reweighted coef- 
ficients, which replace the ill-determined coeffi- 
cients, thus converge to accurate values. 

3.2. Results 

Measurements of the chiral (ipip) and diquark 
condensates (ipip) at strong ]T^] and intermediate 
coupling p7[|18|[19[|[p0[, indicate at least two dis- 
tinct finite density regimes. The first, a low den- 
sity vacuum regime in which (ipip) ^ (ipip) = 0, 
is believed to be separated from a thermodynamic 
regime, where (ipip) = (ipip) ^ 0, by a transition 
at ll — im^. The U(2) symmetry at m, ll = al- 
lows a rotation between the two condensates and 
so the superfluid diquark can be used as an indica- 
tor of the spontaneous breaking of the staggered 
fermion chiral symmetry remnant U(1)a (in ad- 
dition to U(l)v). 

Our measurements of the zeros of the Grand 
Canonical Partition function identify this transi- 
tion as we tabulate in Table f . By evaluating the 
Lee Yang zeros a n in the complex /i plane we can 
readily identify the value of ll at the transition. 
In the thermodynamic limit \mr\\ — > and so 
Herji corresponds to ll c . 



i-j n = T\na n 



(15) 



Before our composite reweighting procedure is 
applied, 771 is inconsistent and strongly depen- 
dent on the value of ll d used in the reweight- 
ing. However the jacknife error estimates with 
our new scheme indicate that the zeros can be 
consistently determined where we alleviate the in- 
accuracies caused by reweighting. As we increase 
the number of ensembles generated at succesive 
values of ll used to cover the expansion, Re rji 
becomes consistent. Also Im^i — > 0, from which 
we are able to more confidently associate Re rji 
with /i c . Our measurements at j3 = 1.5, m — 0.05 
on a 6 3 4 lattice indicate a transition at [i = 0.48, 
which is in approximate agreement with existing 
measurements fl| || [|| §| || . The prominent 
peaking of our further measurements of the quark 




Figure 1. Lee Yang zeros evaluated in the com- 
plex \x plane for SU(2) at intermediate coupling 
from seven composite reweighted ensembles. 
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Figure 2. ( Re (e 1 ^ 71 )) \ Cn \ the ensemble-averaged 
real part of the phase of the polynomial expansion 
coefficient c„ (for n = 3), which we use as a Monte 
Carlo measure for static SU(3) at intermediate 
coupling. 



number density susceptibility from the composite 
weighted expansion coefficients, we believe, indi- 
cates that this transition is first order. 

4. Static SU(3) at Intermediate Coupling 

Having developed our new approach in SU(2), 
we now evaluate the ensemble-averaged polyno- 
mial expansion coefficients of static SU(3) |25j, 
from canonical ensembles generated with differ- 
ent numbers of background quark sources. In this 
way we are able to again systematically vary the 
localised region in which the ensemble-averaging 
is reliable, and combine coefficients from different 
ensembles. 
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Our implementation of static SU(3) involves 
setting G = in Eqn.(H) which differs slightly 
from existing schemes [£5| ||(| [^7| in that we ad- 
ditionally incorporate the relativistic effect of an- 
tiquarks. The motivation for choosing this model 
over full SU(3) is purely the computational expe- 
diency that refreshing c n in one dimension alone 
affords the Monte Carlo sampling. Although the 
expansion coefficients are in general complex, our 
measurements of the ensemble-averaged real part 
of the phase of our measure indicate that the sign 
problem of the measure is negligible, Fig. (2). As 
n moves away from the index of the expansion co- 
efficient used as the Monte Carlo measure, how- 
ever, the ensemble averaging of the expansion co- 
efficients becomes less effective, as we monitor in 
Fig. (3). For the two ensembles generated with the 
measures C75 and C51, W 75 51 (n) is only reliably 
determined for 50 < n < 65. Elsewhere error es- 
timates indicate that the measurement is unreli- 
able, and so only the expansion coefficients in this 
small range are effectively sampled. Our measure- 
ments at j3 = 5.71, to = 0.1 with a 6 3 4 lattice vol- 
ume of the zeros Fig. (4), and corresponding quark 



number density susceptibility, indicate a first or- 
der finite density transition at /i = 0.098(1). 

5. Conclusions 

Although we have addressed the pathologies 
of the Glasgow method with our new compos- 
ite reweighting approach, the main stumbling 
block to the evaluation of finite density SU(3) 
with dynamical fermions would appear to remain 
the numerical effort required to evaluate the 
reweighting method. We seem to have traded the 
overlap problem for the efficiency our evaluation 
of canonical ensembles, which we believe now re- 
quires parallel computing. However, rootfinders 
can often find the zeros of a polynomial from 
the first few coefficients of an expansion, which 
we can achieve in our method by implementing 
a shift variable to truncate the polynomial pi} ). 
Further progress in evaluating SU(3) with dy- 
namical quarks at finite density can therefore 
be made by evaluating ensembles weighted with 
the first few expansion coefficients alone in this 
manner. 

Thanks to M. Alford, and I. M. Barbour for 
useful discussions. 
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Figure 4. Lee Yang zeros evaluated in the com- 
plex [x plane from eleven composite weighted en- 
sembles for static SU(3) at intermediate coupling. 



